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"q i Abstract: It is expected that the magnetic field in the surface of magnetars do not exceed 

2 ! 10^'^ G. However, in heavy ion cohisions, this value is expected to be much higher. We 

investigate the effects of a magnetic fields varying from 10^^, to 10^*^ G in strange matter 

^\i ' (composed of u, d and s quarks). We model matter as a free gas of baryons and mesons 

K*" ' under the influence of an external magnetic field. We study the effects of such strong fields 

j~^ ■ through a x fit to some data sets of the STAR experiment. For this purpose we solve 

the Dirac, Rarita-Schwinger, Klein-Gordon and Proca equations subject to magnetic fields 

in order to obtain the energy expressions and the degeneracy for spin 1/2, spin 3/2, spin 
r~^ \ and spin 1 particles, respectively. Our results show that a field of the order of 10^^ G 

CN ' produces an improved fitting to the experimental data as compared to the calculations 

without magnetic field. 



Contents 

1 Introduction 1 

2 Formalism 2 

2.1 Spin 1/2 Baryons 3 

2.2 Spin 3/2 Baryons 4 

2.3 Spin Mesons 4 

2.4 Spin 1 Mesons 5 

2.5 Thermodynamics 5 

2.6 Chemical Potential 7 

3 Results and Discussions 7 



1 Introduction 

According to Quantum Chromodynamics, the quark-gluon plasma (QGP) phase refers 
to matter where quarks and gluons are believed to be deconfined and it probably takes 
place at temperatures of the order of 150 to 170 MeV. In large colliders around the world 
(RHIC/BNL, ALICE/CERN, GSI, etc), physicists are trying to find a QGP signature 
looking at non-central heavy ion collisions. 

Possible experiments towards this search are Au-Au collisions at RHIC/BNL and Pb- 
Pb collisions at SPS/CERN, where the hadron abundances and particle ratios are used in 
order to determine the temperature and baryonic chemical potential of the possibly present 
hadronic matter-QGP phase transition. 

In previous papers a statistical model under chemical equilibration was used to calcu- 
late particle yields [1, 2] and in these works the densities of particles were obtained from 
free Fermi and Boson gas approximations, where the interaction among the baryons and 
mesons were neglected. More recently, relativistic nuclear models have been tested in the 
high temperature regime produced in these heavy ion collisions. In [3, 4] different versions 
of Walecka-type relativistic models [5] were used to calculate the Au-Au collision particle 
yields at RHIC/BNL and in [6] the quark-meson-coupling model [7-9] was used to calcu- 
late this reaction results and also Pb-Pb collision particle rations at SPS/CERN. In all 
cases 18 baryons, pions, kaons, p's and K*s were incorporated in the calculations and a 
fit based on the minimum value of the quadratic deviation was implemented in order to 
obtain the temperature and chemical potential for each model, according to a prescription 
given in [1]. For Au-Au collision (RHIC) these numbers lie in the range 132 < T < 169 
MeV and 30.5 < fiB < 62.8 MeV and for Pb-Pb collision (SPS), 99 < T < 156.1 MeV and 
167.5 < HB< 411 MeV. 



On the other hand, the magnetic fields involved in heavy-ion collisions [10-12] can 
reach intensities even higher than the ones considered in magnetars [13, 14]. As suggested 
in [10-12] and [15-17] it is interesting to investigate fields of the order of eB = 5 — 30m^ 
(corresponding to 1.7 x 10^^ — 10^*^ Gauss) and temperatures varying from T = 120 — 200 
MeV related to heavy ion collisions. In fact, the densities related to the chemical potentials 
obtained within the relativistic models framework, in all cases, are very low (of the order 
of 10~^ fm~'^). At these densities the nuclear interactions are indeed very small and this 
fact made us reconsider the possibility of free Fermi and Boson gases, but now under the 
infiuence of strong magnetic fields. 

In a recent paper [18], the author studies the synchrotron radiation of gluons by fast 
quarks in strong magnetic fields produced in heavy ion collisions and shows that a strong 
polarization of quarks and leptons with respect to the direction of the magnetic field is 
expected. The polarization of quarks seems to be washed out during the fragmentation 
but this is not the case of the leptons. The observation of lepton polarization asymmetry 
could be a proof of the existence of the magnetic field, which may last for 1 — 2 fm/c. This 
slowly varying magnetic field could leave its signature in the particle yields. 

The purpose of the analysis we present in this paper is to check if the inclusion of 
strong magnetic fields can improve the fitting of experimental results. We start from the 
simplest possible calculation, assuming that the magnetic field is homogeneous, constant 
and time-independent. We are aware that it is not the case, as shown in [19, 20], where 
the shape of the magnetic field presents a special non-trivial pattern. Moreover, from 
the calculations performed in these references, one can see that after averaging over many 
events one is left with just of the components of the magnetic field. Nevertheless, the event- 
by-event fluctuation of the position of charged particles can induce another component of 
the magnetic fleld (perpendicular to the remaining one in the average calculation) and 
also an electric fleld, which is quite strong at low impact parameters. While the magnetic 
field remains quite high in peripheral collisions, the opposite happens with the electric 
field. To make our first analysis as simple as possible, we shall restrict ourselves to data at 
centralities of the order of 80%, i.e., high values of the impact parameter 6 ~ 11 — 13 fm, 
where we are more comfortable to disregard the electric field effects. 

In the present paper we briefly revisit the formalism necessary for the calculation of 
particle densities subject to magnetic fields and the expressions used to implement a x^ fit 
to the experimental results. 



2 Formalism 

We model matter as a free gas of baryons and mesons under the infiuence of a constant 
magnetic field. We consider only normal and strange matter, i.e., the baryons and mesons 
constituted by u, d and s quarks: the baryon octet (spin 1/2 baryons), the baryon de- 
cuplet (spin 3/2 baryons), the pseudoscalar meson nonet (spin mesons) and the vector 
meson nonet (spin 1 mesons), which leaves us with a total of 54 particles (18 baryons, 18 
antibaryons and 18 mesons). 



We utilize natural units {h = c = 1) and define eo = /xq = 1. From the relation 
we obtain that the electron charge is e = \/47ra, where a = -r^ is the fine 



structure constant. The natural units with the electron charge in that form is known as 
Heaviside-Lorentz units [21]. 

In this work, the magnetic field is introduced trough minimal coupling, so the deriva- 
tives become 

df,^D^ = d^+iqA^. (2.1) 

We write the charge as q = eq\q\, where €q = +(— ) corresponds to a particle with 
positive (negative) charge, and assume the gauge 

A^ = 6^2XiB -^ Ao = and A = {0,xiB,0), (2.2) 

so, 

V-A = and VxA = Be3, (2.3) 

and the derivatives 

Df, = df,- ieq\q\Bxi6f,2. (2.4) 

We search for solutions of the fields ip in the form 

(.)^f d^)e-^*+-^^-^ (, = 0) 

where ipa are the components of the field ip and e = +(— ) corresponds to the states of 
positive (negative) energy. 

For the spin 1/2 baryons (Dirac field) ip has 4 components, for the spin 3/2 baryons 
(Rarita-Schwinger field) ip^ has 16 components, for the spin mesons (Klein-Gordon field) 
ip has just one component, and for the spin 1 mesons (Proca field) V^ has 4 components. 

Due to the use of statistical methods to deal with the system under consideration, we 
do not need the complete expression for ip, but just the form of the energy E for each one 
of the fields and the degeneracy of the energy levels 7. 

2.1 Spin 1/2 Baryons 

The baryons with spin 1/2 are described by the Dirac Lagrangian density [22] 

CD = ^[iYD^ - m)^, (2.6) 

which (after we apply the Euler-Lagrange equation) lead us to the equation of motion 

(i7^L»^ - m)tp = 0. (2.7) 

where 7^^ are the Dirac matrices. 

The solution of the equation of motion gives 



^ ^ f \/p^ + 'rm? {q = 0) ,^ 

\ ^/pi + m? + 2v\q\B (g / 0) ' ^ ' 



where v runs over the possible Landau Levels and the degeneracy for the energy states are 
given by: 

7=1 ' ^' = °^ (29) 

2.2 Spin 3/2 Baryons 

The baryons with spin 3/2 are described by the Rarita-Schwinger Lagrangian density [23, 

24] 

C'''' = -^V5/.(e^"''"757.^p + iman^^, (2-10) 

where 7^ = i'^'^^^'^'^^^ and cr^'^ = §[7^,7'^]. 
The equation of motion reads 

{i-f^'D^ - ?n)V'. = with 7^^^ = and L^'Vm = 0- (2-11) 

The solution of the Rarita-Schwinger equation is not trivial and poses non-causality 
problems. To obtain the degeneracy of the energy states, we follow the prescription used 
in [22], which is given in detail for the Rarita-Schwinger equation in [25]. Observing the 
equation of motion one can see that each component of ip^ obeys a Dirac type equation, 
so the energy must have the form 



E=\ ,,^^t""\ , i^ = °?. (2.12) 

Besides that, V'^ has 4 components, but, two equations are constrainted, which means 
that only 2 components of 'i/'/j are really independent. So, V'/x have 2 polarizations, but, 
(because of the Dirac equation solution) each polarization are double degenerate. In the 
presence of a magnetic field there is another constraint for the 1/ = and v = \ energy 
levels, which leads to the following degeneracy for the energy states 

^=\ 4-25,o-<5.i (g/O)- ^^-^^^ 

2.3 Spin Mesons 

The mesons with spin are described by the Klein-Gordon Lagrangian density [26] 

r^c = D^^*D^^ _ mV>, (2.14) 

whose equation of motion is given by 

(D^'Z?^ + m^)V' = 0, (2.15) 

with the energy satisfying the relation: 



■sjfp^+rn?' {q = 0) 

^pi+m? + {2u + l)\q\B {q / 0) 



(2.16) 



2.4 Spin 1 Mesons 

The mesons with spin 1 are described by the Proca Lagrangian density [27] 

The equation of motion is 

{D^'D^ + m^)il)y = with D^ip^" = 0. 



(2.17) 



(2.18) 



Each component of ipn obey a Klein-Gordon type equation, so that the energy states 



are 



E 



VW+ 



(9 = 0) 

^pi + m? + {2u + l)\q\B {q ^ 0) ' 



(2.19) 



ipfj^ has 4 components, but, one of the equations is a compressed constraint equation, 
which means that only 3 components of ip^ are independent. So, each energy state have 
3 polarizations in the case with zero charge (or without magnetic field). If the charge is 
different from zero (and we have the presence of an external magnetic field) there is an 
additional constraint for the v = energy level, which leads to the following degeneracy 
for the energy states 

f 3 iq = 0) 
3 - 5,0 (g / 0) ■ 
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(2.20) 



2.5 Thermodynamics 

Using the Grand Canonical formalism we obtain that the particle densities for the baryons 



are 



Pb= { 



1 /'°° 
7fe^ / ■^'^^^ ~ l-ib)p^dp {q = 0) 

E^^^ / fiEb-^^b)dp (g/o) 

i/=0 



2^2 7o 



(2.21) 



for the antibaryons are 



1 f°° 
^b^ J f{Eb + fib)p^dp iq = 0) 



oo 
v=0 



qb\B '■°° 
27r2 



fiEh + fib)dp (q^O) 



(2.22) 



and for the mesons are 



1 /'°° 
7m ^ / ^(^"^ ~ ^im)p^dp {q = 0) 



u=0 



\qm\B 



(2.23) 



m a / ^(^"» ~ t^m)dp {q / 0) 



with f{x) = (e^/^ + l)"i and h{x) = (e^/^ - l)-^ . 



The total baryonic particle density is 



PB ='^iPb- Pab) 
b 



(2.24) 



and the total mesonic density is 



PM = ^ Pm- (2.25) 

m 

The energy density is given by the sum of the energy densities of each particle, so 



e = '^{eb + f^ab) + X^e^ 



(2.26) 



with 



^b= < 



1 



767TZ? / Ebf{Eb - Hb)p^dp {q = 0) 
qb\B 



oo 
i/=0 



poo 

/ Ebf{Eb-iib)dp (g/0) 
Jo 



1 f°^ 
7b^y EbfiEb + fib)p^dp iq = 0) 

°° |„ I D /-OO 

Y^ 7b^ / EbfiEb + fib)dp {q / 0) 



1 

'2^ 



5^ 7m 



27r2 



Emb{E„^ - llrri)P^dp {q = 0) 
/■oo 

/ Emb{Em - Pm)dp {q / 0) 
JO 



in the same way the pressure is given by 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



with 



Pb= { 



1 f°° 1 
7fe^ / -^f{Eb- Pb)p^dp iq = 0) 



E 

i/=0 



76- 



27r2 



1 



fiEb-fMb)p^dp (g/0) 



P, 



ab 



1 1"°° 1 

76^ y —f{Eb + Pb)p^dp (g = 0) 

Y^^b-fr YfiEb + Pb)p'dp (g/0) 

v=0 '^ -JO b 



Prr 



1 



1 



oo 
iy=0 



7m7r3 / -^h{Ein - Hm)p^dp {q = 0) 

JO ^m 

\Qm\B 

27r2 



/■oo -|^ 

/ -=^b{Em - Pm)p'^dp {q / 0) 

JO ^m 



(2.31) 



(2.32) 



(2.33) 



the entropy density can be found through 

S = e + P -"^flbipb- Pab) -^f^mPm- (2.34) 



m 



2.6 Chemical Potential 

The hadron chemical potential is 

Ph = Bh PB + hh Ph + Sh PS, (2.35) 

where Bh, /3/j and Sh, are respectively the baryonic number, the third isospin component 
and the strangeness of the particle h. The baryonic chemical potential pB is a free param- 
eter of the system (the other is the temperature T) . The chemical potential of isospin pj^ 
and strangeness ps are determined trough their respectively conservation laws. 

We impose the local conservation of the baryonic number, isospin and strangeness. 
This imposition leads to the following equations 

^Bhph = -^, ^hhPh = y, ^Shph=y, (2.36) 

h h h 

where Nb is the total baryonic number, I^ is the total isospin, S is the total strangeness 
of the system and V are the volume occupied by the system. The charge conservation is 
automatically achieved trough the other three conservation laws. 

The baryonic number of an Au atom is Nb = {N + Z) = 79 + 118 = 197, the isospin is 
/3 = (Z - N)/2 = 19.5 and for the deuteron (d) we have that A^^ = 1 + 1 = 2 and I3 = 0. 
Hence, assuming that the total strangeness of the system is zero, we write the following 
table for the conserved quantities: 

• Au+Au Collision, Nb = 394, I3 = -39, 5 = 0. 

• d+Au Collision, Nb = 199, I3 = -19.5, 5 = 0. 

At this point it is important to emphasize some of the drawbacks of our simple calcu- 
lation. As shown in [28], the magnetic field should depend on the charges of the colliding 
nuclei and the number of participants should vary for different centralities. These con- 
straints were not taken into account directly in our calculations. All the information we 
use as input come from the experimental particle yields and the magnetic field is modified 
until the best fitting is encountered. The number of different participants is reflected only 
in the resulting radii. 

3 Results and Discussions 

We have implemented a x^ fit in order to obtain the temperature and chemical potential. 
The particle properties (spin, mass, baryonic number, isospin and strangeness) were taken 
from the Particle Data Group [29]. 



In tables 1, 2, 3 and 4 we show our results corresponding to the temperature and 
chemical potential that give the minimum value for the quadratic deviation x^: 

f^exp q~itheo\2 
i * 

where 7^^^^ and TZf^^° are the i particle ratio given experimentally and theoretically, and 
(Tj represents the errors in the experimental data points. 

To make clear the improvement in the data fitting by the addition of the magnetic 
field, we calculate the relative percent deviation (A%) with respect to the experimental 
values for B = and the best B ^ (the bold columns in the tables) trough the equation 



A 



% 



'jytheo jyexp 



'jyexp 



100%, (3.2) 



and show these values in parenthesis in all the tables. 

For the simulations our code deals with 5 unknowns (/^s, /^/s, /-f^, T, V) and 3 con- 
strained equations. We run over the values of ^ib and T (the free parameters) in order to 
find the smallest x^- Our results are given next. 

In tables 1, 2, 3 and 4, B is the magnetic field, T is the temperature, ^b is the baryonic 
chemical potential, x^ is quadratic deviation, /i/3 is the isospin chemical potential, ^s is 
the strangeness chemical potential, R is the radius of the "fire-ball", pB = J2biPb ~ Pah) is 
the usual baryonic density, p^ = p^++ — Pa++~^Pa+ ~ Pa+ +Pao ~ Pa" +Pa- ~ Pa- i^ delta 
baryon density, pM = X^^ Pm is the meson density, p^ = p^o +Pt^+ +Pn- is the pion density, 
e is the energy density, P is the pressure, s is the entropy density and ndf is the number of 
degrees of freedom. For i? = 0, ndf = 5 (7 experimental values minus 2 free parameters, T 
and /i), for S 7^ 0, ndf = 4 (7 experimental values minus 3 free parameters, T, p, and B). 
7r~/7r^, K~ /K^, p/p, K^ /t:~ , K'^ jix^ and p/vr"*" are the theoretical (first 7 columns) and 
experimental (last column) particle ratios [30]. The temperatures and baryonic chemical 
potentials obtained from the statistical model in [30] are also given in the last columns of 
all tables. 

In figs, l-a/6, 2-a/h, 2>-a/h and A-a/h we plot the experimental and theoretical ratios 
for B = {) and the best i? 7^ 0. In figs. 1-c, 2-c, 3-c and 4-c we show the y^ behavior for 
B = Q and for the best B ^ 0. In figs. 1-d, 2-d, 3-d and 4-d we show the x^ behavior for 
the different magnetic fields. One can notice that the best fitting is generally obtained for 
magnetic fields around 6 m^, a little higher than what is expected for RHIC collisions (5 



ml] 



Our results show that, even for the free Fermi and Boson gas models, a strong magnetic 
field plays an important role. The inclusion of the magnetic field improves the data fit up 
to a field of the order of B = 10^^ G. For stronger magnetic fields, it becomes worse again. 
This behavior is easily observed in tables 1 to 4 and in figs. 1-d to 4-d. It is worth pointing 
out how the "fireball" radius R and the total density p vary with the magnetic field in a 
systematic way: R and p practically do not change between B = and B = 10^® G, but 
when the field increases even further, the density increases and the radius decreases. This 



behavior is common to all collision cases studied. This huge jump in the density explains 
why the ratios get worse for a magnetic field of the order of B = lO^'^ G, for which the 
densities are much higher than what is expected in a heavy ion collision. 

Our model gives a good description for the particle/antiparticle ratios, but fails to 
describe the relation between baryons and mesons. This occurs because our model produces 
too many mesons (especially pions) as shown explicitly in the particle densities. In all 
collision types our model presents a baryon density (pb) with more than 30% of A baryons 
and a meson density (pm) with more than 60% of vr (pn)- The relative percent deviations in 
the particle yields show clearly that some results improve considerably when the magnetic 
field is considered, while others remain unaltered or even get slightly worse. However, 
our figures also show that the behavior of the x^ changes drastically with the addition of 
the magnetic field and that the temperature and chemical potentials calculated with the 
statistical model lie within the 3 — a confidence ellipse obtained for the best x^ in some 
cases, but they are always outside the confidence ellipses obtained with zero magnetic field. 

We would like to comment that when we first started these calculations, we were not 
aware of references [19, 20] and we used data obtained for low centralities, i.e., low impact 
parameters. In that case, the minimum x^ was generally smaller than the ones shown in 
this work and we believe this was so because of the larger error bars accompanying data 
at low centralities. 

Further improvements on the presented calculations are under investigation, namely, 
the inclusion of electric fields at low impact parameters and the variation of both electric 
and magnetic fields with the number of participants in the collisions. Moreover, we are 
working on the inclusion of the anomalous magnetic moments and in the description of 
pion-pion interactions. We next intend to repeat these calculations for the ALICE/LHC 
data for the future Au+Au runs with all these improvements, so that our results become 
more realistic. 
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Table 1. Au+Au (70-80%) -y/sjvw ~ 200 GcV. Results obtained for different values of the magnetic field. 
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Figure 1. Au+Au (70-80%) collision at y^ijviv = 200 GeV. (a) particle/ antiparticle ratios. (6) 
mixed ratios, (c) and {d) x^ behavior. 
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Table 2. Au+Au (58-85%) collision at y/s^N ~ 130 GcV. Results obtained for different values of the magnetic field. 
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Figure 2. Au+Au (58-85%) collision at y/TNN = 130 GeV. (a) particle/ antiparticle ratios. (6) 
mixed ratios, (c) and (d) ^ behavior. 
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Table 3. Au+Au (70-80%) collision at -y/sjvA?" ~ 62.4 GcV. Results obtained for different values of the magnetie field. 
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Figure 3. Au+Au (70-80%) collision at ^^snn = 62.4 GeV. (a) particle/antiparticle ratios. (&) 
mixed ratios, (c) and (d) x^ behavior. 
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Table 4. d+Au (40-100%) collision at -y/sjvw = 200 GcV. Results obtained for different values of the magnetie field. 
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Figure 4. d+Au (40-100%) collision at y^ijvTV = 200 GcV. (a) particle/antiparticle ratios. (6) 
mixed ratios, (c) and (d) x^ behavior. 
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